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We introduce the concept of the a point of minimal capacity of the domain, and observe 
a connection between this point and the lowest eigenfunction of a Laplacian on a domain, 
in one special case. 

1 A physical motivation and precise definitions 

1.1 The "warmest" point of a domain. 

For a given a domain in D C M. n , let m(D) G D be a point of maximum of the lowest 
eigenfunction of the Laplacian with zero Dirichlet boundary condition^. A point m(D) is 
in a certain sense the warmest point in D: if the domain, viewed as a heat-conducting 
medium, starts with a positive temperature distribution and is cooled by the maintenance 
of zero temperature on the boundary, then after a long time a (local) maximum of the 
temperature will approach a maximum of the lower eigenfunction. Indeed, the solution u 
of the heat equation 

u t = An, u(0,x) = u (x), u(t,x)\ xedD = 

on a domain D in 1" is given in terms of eigenfunctions v n of the Laplacian with the 
Dirichlet boundary conditions: 



,t) = ^2a n e Xnt v n (a 



u[x, 

n=0 



with > Ao > Ai > • • • . The leading mode u(x, t) = e A °*t>o(x) becomes dominant for t large 
(we assume ao > and normalize to ao = 1). Hence the warmest point indeed approaches 
a maximum of uq, as claimed. For the background on the Dirichlet problem we refer to the 
classical text [I]. 



1.2 The minimal capacity point in 2D. 

Along with the maximizer m(D) of the lowest eigenvalue we consider another special point 
which we will call the minimal capacity point of the domain D. In a sense made precise 
shortly, this is the "best insulated from the boundary" point; we define it as follows. Let 
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1 Such a point need not be unique for non-convex domains. 
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Figure 1: The minimal capacity point m minimizes the flux from the small disk around m 
to dD; a more precise definition is in the text. 

S e (x) C D be the sphere of radius e centered at x & D, Figure [TJ Let us prescribe the 
boundary conditions 

«S £ (x) = 1, and u 9D = 0, 

and let u be the solution of the resulting Dirichlet problem on the cored domain D \ 
mt(S e (x)). If we interpret u as the temperature, heat will flow from the warm sphere S £ (x) 
to the cold boundary dD, with the heat flux given by 

F E (x) = - [ Vu-ndS. (1) 

This flux depends on the choice of x and on e. Let x e be a minimizer of F e (such a minimizer 
may not be unique, for instance for a dumbbell-shaped domain with a thin neck). If 

c(D) = f lim x £ 

exists, we will call this limit the minimal capacity poin^. Intuitively, c(D) is the point best 
insulated from the boundary, since it minimizes the flux of heat. 

An alternative physical interpretation: c{D) is the point of least capacity of the capacitor 
whose electrodes are dD and S x {e), in the limit of e — > 0. Indeed, if we interpret Vu as the 
electrostatic field in the vacuum (which we can do since Au = 0), then the flux F £ (x) of this 
field is, by Gauss's law, and up to a scaling factor, the amount q of electrostatic charge on 
the spherical electrode S x (e). But the potential difference V between the electrodes S e (x) 
and dD is V = 1, and by the definition of the capacitance (C = q/V = F £ (x)/1 = F £ (x)) 
we conclude that F e (x) is precisely the capacity of the capacitor in question. 

There is yet one more interpretation of F £ {x) defined by Eq. ([1]): it is the potential 
energy (up to a factor of |) stored in an elastic membrane subject to the boundary conditions 
mentioned in the definition of the least capacity point (for small displacements for which 
the nonlinear equation of the minimal surface can be replaced by Vu = 0). One can imagine 

2 No claim is made as to uniqueness of c(D); in fact, one can construct examples with multiple such 
points. In the case we consider, however, c is unique. 
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placing a heavy slippery disk S £ on the horizontal membrane. The disk will slide to the 
position of least potential energy ([I]), i.e. to the vicinity of the least capacity point. 

To summarize this and the preceding sections, we mentioned two physically reasonable 
definitions of the "warmest point" in D. The vague intuitive connection between these two 
definitions suggests a more precise mathematical relationship. The goal of this note is to 
explore this relationship on a simple example (Section [2]), and to state an open problem 
(Section [3]). 



We mention in this connection another problem suggested by Walter Craig and studied 
in depth by Jochen Denzler [2] [3]. This problem addresses the question "where to place 
a window of given area to minimize heat loss?" The problem reduces to minimizing the 
principal eigenvalue of a mixed Neumann-Dirichlet problem. 

1.3 A characterization of the least capacity point. 

Before discussing the main result we give yet one more characterization of the minimal 
capacity point c(D) in terms of Green's function of the domain. 

Theorem 1. Let D C W 1 be a domain for which any Dirichlet problem has a solution. 
Given any p E D let v p (z) be the harmonic function of z satisfying Dirichlet boundary 
conditions v„ = — In 1 on dD if n = 2 and v„ = —-, — t^tt on dD if n > 2. The least 

P \z—p\ J P \z—p\ n 1 J 

capacity point c(D) maximizes the function v z (z). 

Remark. The definition of v p {z) is motivated by the desire to make Green's function 

vJz) + In - — - — - 
\z — p\ 

(for n = 2) vanish on dD. 



Proof of Theorem 1. We wish to construct the harmonic function involved in the defini- 
tion of c(D). As a candidate, we take Green's functiorH 

u(z) = k- l ( \n\z -p\- 1 + v p (z) ), (2) 



where p G D and where k will be chosen so that the average 

U\z~p\=e = 1- (3) 

By the definition of v p we have uqd = 0. 

Averaging Eq. ([2]) over the sphere S e (p) and using the fact that v p (z) is harmonic, we 
obtain the condition on k which guarantees ([3]): 

l = k- 1 Qnl/s + v p (p)), 



3 we consider the case n = 2, leaving out the obvious changes required for n > 2. 
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or 



k = {\nl/e + v p {p))- 1 . 
Thus the harmonic function satisfying uqd = and U| 2 _ p | =e = 1 is given by 

_ In \z-p\~ 1 + z; p (z) 
HZ) ~ lnl/e + v p (p) ' 

to the leading order for small s. Computing the "heat flux" we get the contribution 2tt 
from the logarithmic term and zero from v p {z) since the latter is a harmonic function in D; 
to the leading order we have 

/ 2-7T 

F £ (p) = — <f> Vn • nds 



J \z— p\=e 

lnl/e + v p (p)' 

This shows that for small e, the minimizer x £ of F £ (p) is close to the maximizer of v„(p). 
In other words, c(D) = lim^o^e is the maximizer v p (p). 

1.4 The point of minimal capacity for an elliptic operator on an interval. 

For the case of n = 1 the concepts described before become trivial: both points m{D) and 
c(D) are simply the midpoints of the interval D. However, the question is still interesting 
for a more general elliptic operator 

Lu = (a(x)u)' , a{x) > (4) 

on C 2 [0, 1] with the homogeneous boundary conditions u(0) = u(l) = 0. 

Since this is not a special case of the above, we repeat the definitions of the points m 
and c. Since the domain can be taken as the unit interval [0, 1], the points will essentially 
depend only on the operator L, and we will write c[L] and m[L] in slight break with the 
earlier notation. 



Definition of c[L]. Let us introduce Green's function: for any < s < 1 (Figure [2]) 
we consider the solution u-(x;s) of Lu = on x € [0,s] with u-(0;s) = 0, n_(s;s) = 1. 
Similarly we define u+ as the solution of Lu = on [s, 1] with m + (s;s) = 1, u+(l;s) = 0. 
We consider the (one-dimensional) flux out of x = s, the analog of ([1]): 

F(s) = a(s)(-u' + (x; s) + u'_(x; s)) x=s , (5) 

where ' denotes the x-derivative. Note that the signs in ([5]) are chosen so that the flux out 
of x is counted with a positive sign, just like in the higher dimensional case ([T]). Finally, 
we define c = c[L] as the minimizer of F(s). It is clear that such a minimizer exists, since 
F(s) — y oo as s — ^ or s — t 1. 

Remark 1. The least capacity point c = c[L] is alternatively defined by 

a^ 1 (x) dx = J a^ 1 (x) dx. (6) 
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Figure 2: F(s) > F(c): the flux from c is smallest, i.e. c is the best insulated point. As 
it turns out, c also bisects electrical resistance of the segment if a(x) is interpreted as the 
local conductivity. 



In particular, the least capacity point is unique. 



Proof. From the definition of u± as the solutions of Lu = with appropriate boundary 
conditions we obtain 

Jo a Hi)** ^a- x {t)dt 

Inserting these into © we obtain 

1 1 1 1 

F s ) = — — + 



/ a-i(t)dt f;a-Ht)dt R(s) R(l)-R(s)' 

where = J S a~ 1 (t) dt. Now since f(r) = + is maximized by r = k/2, we conclude 
that F(s) is maximized by that value of s which gives R(s) = i.e. by s = c satisfying 

®. " 



If a(x) is interpreted as the local electrical conductivity of a wire at x (that is, the 
conductivity per unit length at x), then R(s) = / Q s a~ 1 (t) dt is the resistance of the piece 
of wire [0, s]. Note that (jSJ) is intuitively plausible: it states that for c to maximize the 
electrical resistance from itself to the two ends of the segment, c must bisect the resistance 
of the entire segment. One can also give an equivalent thermal interpretation of (jBJ) by 
replacing the word "voltage" by "temperature", "current" by "heat flux", etc. A mechanical 
interpretation of (jSJ) is similarly simple: it states that c is the point on a string for which the 
pieces [0, c] and [c, 1] have equal Hooke's constants. By the definition, on the other hand, c 
is the point which is "easiest" to slide in the x-direction if a spring is grabbed at c. 



Remark 2. Operator @ arises in many physical settings of which we describe briefly 
three. 

1. Transversal vibrations of a string with variable linear density p(x) are governed (to 
the leading order) by the wave equation 

p(x)u u = u xx . 
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By introducing the mass parameter s via s = Jq p(y) dy we rewrite the above ODE 
in the form 

d 

u t t = (a(s)u) , a{s) = p(x(s)), ' = —. 

2. Longitudinal vibrations of an inhomogeneous string. Consider an elastic string fixed 
at two ends and undergoing (small) longitudinal vibrations along the x-axis, with the 
ends of the string fixed. Let u(x, t) denote the displacement from the equilibrium of 
that particle of the string whose equilibrium position is x. Assuming that the material 
of the string satisfies linear stress-strain relationship, we can rewrite Newton's second 
law as 

^tEMY, •=», ' = (8) 

where E is Young's modulus, p = p(x) is the linear density of the string and A = A(x) 
is the variable cross-sectional area of the string. By choosing the new space variable 
s just as in the preceding example, one reduces (|8|) to the same form as above: 

u = (au')'. (9) 

One can interpret the last equation directly as describing longitudinal vibrations of a 
string with constant linear density p = 1 and with variable "local" Hooke's constant 
a = a(x). 

3. The parabolic PDE 

u = {au')' (10) 

is the heat equation describing the evolution of temperature u along a rod (with 
insulated walls) with heat conductivity a(x) and with heat capacity (specific heat) 
equal to one unit per unit lengtlJll. Indeed, interpreting a as the heat conductivity 
amounts to saying that —au' is the heat flux along the rod at x. The instantaneous 
rate of heat gain by a segment [x, x + dx] is then (-au')'dx, to the leading order; since 
the heat capacity is 1, this results in (fTUj) . 



4. The ODE 

Lu = (au)' = 

governs an electrostatic potential u along a resistive wire with local conductivity a{x). 

x+dx 

Indeed, — a(x)u(x)' is the current through point x, so that (au') = expresses 

X 

the conservation of charge in a segment [x,x + dx] (Kirchhoff's first law). 



2 The main result. 

In the following theorem we consider the operator L given by (|4|). As before, we denote by 
c = c[L] the point of least capacity and by m = m[L] the point of maximum of the principal 
eigenfunction ("the warmest point"). 

4 one can reduce the equation with a variable heat capacity to this one by the same transformation as in 
the preceding example. 
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Theorem 2. Consider the elliptic operator Lu = (an')' on [0, 1], where a > is a monotone 
increasing -function on the interval. Then c[L] < m[L], where c[L] is the point of least 
capacity associated with L, as defined above, and where x = m[L] gives the maximum to the 
principal eigenfunction of L. 




c[L] m[L 



Figure 3: For a tapered vibrating string, the most elastic point is closer to the thin end 
than the maximum amplitude point of the fundamental mode. 



Corollary 1. Consider vibrations of an inhomogeneous elastic string, either longitudinal 
or transversal, described by u# = (a(x)u')' . Assume that the "stiffness coefficient" a(x) 
increases from left to right, just as in the last theorem. Then the the point of maximal 
amplitude of the fundamental mode lies between the "softest point " c and the thick end of 
the string. The same holds for any interval between two consecutive nodes of any overtone. 

Figure [3] illustrates this effect for transversal oscillations. 



Proof of the theorem. Let u{x) be the first eigenfunction of the operator ([4]): 

(a(x)u'Y = -An, A > 0, (11) 

which means that u > on (0, 1) and satisfies the homogeneous boundary conditions. We 
rewrite this relation system 

[ a(x)u' = v 

, , (12) 
I v = —An 

Let us translate the definitions of m[L] and c[L] into geometrical terms. Since c[L] is defined 
by the condition (jfij), it is natural to choose t = t(x) = Jq a~ 1 (s)ds (the resistance) as the 
new independent variable, so that f)12|) become 

f u = v , 

S • (13) 
1 V = -Xa(x(t))Y, 



where ' = and U(t) = u(x(t)), V(t) = v(x(t)). Assume without the loss of generality that 

a~ 1 (s)ds = 1, 



dt 

■1 







so that t ranges in < t < 1. Then c = c[L] is defined by 
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On the other hand, since x = m is the maximizer of u(x), we have v(m) = u'(m) = 0, i.e. 
V(t(m)) = 0. Introducing the angle 9 = arg(t7 + iV) we restate the definition of x = m as 
9(t(m)) = 0. Since t is a monotone increasing function of x, proving c < m is equivalent to 
proving that t(c) < t(m). Since t(c) = 3, the proof reduces to showing that 

2 < (14) 




0(* m +r*)=0(l) = 7r/2 

Figure 4: Proving that c[L] < m[L], i.e. that 5 < t m . 
Now the angle 9 satisfies 



sin 



2a - \a(x(t)) cos 2 9 = f f(9,t) 



(15) 



along with the boundary conditions 6(0) = 5 and 6(1) = — |. Note that a(x(t)) is a 
monotone increasing function and that / < for all values of its arguments, so that 6 is 
monotone decreasing. 

The key to the proof of (|14[) is the monotonicity of a(x(t)) and the property f(6,t) = 
f(—6,t). The idea is to compare the angle 6 at equal times r before and after crossing the 
zero value, as in Figured! That is, we introduce 



<p(r) = -6(t m - t) and i/j(t) = 6(t m + r), 

so that ip(0) = 4>(0) = 0. Moreover, these 99, ip satisfy 

J ip = 6(t m - r) = f(-(f, t m - t) = f(cp, t m - r) 
1 ip = f(ip,t m +r), 



(16) 



(17) 



where / has been defined in (|15p . Since A(t) = a(x(t) is monotone increasing, f(9,t) is 
monotone decreasing in t and we have f(tp,t m — r) > f(f,t m + r). Since (p, ip share the 
initial condition, the comparison theorem applies: 



7T 



- < il>(r) < <p(r) < 



(18) 
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for r > 0, < t m + t < 1. Thus for some r = t* we have 

-|=^r*)<^(r*)<0 
or, recalling the definition of tp and tp: 

-| = 0(i m + t*)<-0(t m -i*)<O. 

By (^4) we have 

*m + T* = 1. (19) 

On the other hand, from the above we have 9(t m — r*) < 7r/2, which implies t m — t* > 
(since 6*(0) = 7r/2 and 6 < 0). Adding this to (fl~9j) we obtain i m > |, thus completing the 
proof. 

3 An open problem. 

Consider a planar domain D (Figure [S]) bounded by two curves y = ±/(x), where / is a 
positive monotone increasing function, and by two lines x = 0, x = 1. By symmetry, both 
points c(D) and m(D) lie on the x-axis; let us denote by their x-coordinates by c and m 
respectively. 



Open problem 1: show that c < m. 



Open problem 2: show that c < m for the Laplacian with mixed Dirichlet-Neumann bound- 
ary values: Dirichlet on x = 0, x = 1 and Neumann on y = ±/(x). 



A 



y=f(x) 



m 



Figure 5: c < m? 



Remark. If we replace y = f(x) by y = ef(x) with a small s, then the membrane modeled 
by our mixed Dirichlet-Neumann problem becomes so narrow as to resemble an elastic 
string described in Theorem [2j This theorem in fact suggested the second open problem. 
For small e our mixed Dirichlet-Neumann problem is approximated by the one-dimensional 
operator from Theorem [2j There is a vast literature on partial differential operators on thin 
domains; we refer to 0] and to references therein. 
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